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Local continuity equations involving background fields and variations of the fields, are obtained 
for a restricted class of solutions of the Einstein-Maxwell and Einstein- Weyl theories using a new 
approach based on the concept of the adjoint of a differential operator. Such covariant conservation 
laws are generated by means of decoupled equations and their adjoints in such a way that the corre- 
sponding covariantly conserved currents possess some gauge invariant properties and are expressed 
in terms of Debye potentials. These continuity laws lead to both a covariant description of bilinear 
forms on the phase space and the existence of conserved quantities. Differences and similarities with 
other approaches, and extensions of our results are discussed. 
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I. INTRODUCTION 

The study of the conservation laws is present in practically all areas of physics. In the context 
of the theories involving gravity, such a study becomes particularly interesting because of the lack 
of conserved currents representing the conservation of energy and momentum. Moreover, in the 
covariant description of a symplectic structure on the phase space of classical systems, a covariantly 
conserved two-form on such phase space is required. Such descriptions yield a canonical formalism 
manifestly covariant in field theory, which represents a starting point for a quantum mechanical 
treatment . More specifically, in this context, the phase space is defined as the manifold of solutions 
of the classical equations of motion. Such a definition has the advantage of manifest covariance Q . 
The existence of a two-form (say J^), covariantly conserved (9^ — 0) on such phase space, lead 
to a symplectic structure ui defined as uj = J^J'^dJ^^ (where S is an initial value hypersurface) , 
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independent of the choice of S. Moreover, the condition that Poisson brackets satisfy the Jacoby 
identity, holds if u is closed (which means that w is a two- form whose exterior derivative vanishes). 
This condition is satisfied if itself is closed. Several methods have been proposed in the literature 
for constructing these conserved currents, from the guess for the explicit form of the currents Q, 
to the very general result on the existence of a conserved current (the symplectic current) for the 
field variations for any field theory derived from a Lagrangian action ^ . Alternatively, equivalent 
conserved currents have been obtained directly from the equations for the field variations using the 
fact that there exists a conserved current associated with any system of homogeneous linear partial 
differential equations that can be written in terms of a self-adjoint operator ^. Such results 
have been obtained without any statement about the background fields. Nevertheless, the currents 
constructed in these ways are associated with the variations of the gauge fields. For example, in the 
case of the Einstein-Maxwell theory, such currents are expressed in terms of the metric and vector 
potential variations, and in this manner are defined modulo gauge transformations [1-4]. However, 
as it will be shown in the present paper, one can construct, for a restricted class of spacetimes, 
conserved currents (and the geometrical structures associated) without the gauge freedom of the 
matter fields. The approach is based on the concept of the adjoint of a differential operator, which 
represents an extension of that employed in Refs. ^. Our starting point is the existence of 
appropriate decoupled equations for the restricted class of spacetimes under study, which have been 
obtained using the Newman-Penrose formulation jH ^. It is shown that the conserved currents are 
also independent of the perturbed tetrad gauge freedom, which appears when the variations of the 
Newman-Penrose quantities are involved. 

A previous application of the present approach has been given in the study of the perturbations 
of the Schwarzschild black hole, leading up to an important result on the existence of conservations 
laws for non-Hcrmitian systems, which has a direct extension to the scheme of string theory j?]. In 
this manner, another of the purposes of this article is to establish conservation laws for more general 
spacetimes, analogous to that appearing in Ref. 0. However, unlike the results given in that 
reference, one of the main novelties in the present article is the interpretation of the corresponding 
conserved currents as bilinear forms on the phase space, which is connected with the very important 
concept of symplectic structures. This subject will be extended later on. 
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The outline of this paper is as follows. Section II is dedicated to establish the general relationship 
between adjoint operators and conserved currents. The application of this very general result for 
identifying the conserved currents for the field variations in the Einstein-Maxwell and Einstein- Weyl 
theories is given in Sec. III. The properties of such conserved currents, and their relations with the 
existence of symplectic structures on the phase space are discussed. We finish in Sec. IV with some 
concluding remarks and future extensions of the present results. 



II. ADJOINT OPERATORS AND CONSERVED CURRENTS 

New branch: In Rcfs. Q it has been shown that there exists a conserved current associated 
with any system of homogeneous linear partial differential equations that can be written in terms of 
a self-adjoint operator. This result is limited for a self-adjoint system, for which the corresponding 
conserved current depends on a pair of solutions admitted by such a system. However, as we shall 
see below, there exists a more general possibility that extends for systems of equations which are 
not self-adjoint necessarily. The demonstration is very easy, following the basic idea of Refs. 

In accordance with Wald's definition if £ corresponds to a linear partial differential operator 
which maps m-index tensor fields into n-index tensor fields, then, the adjoint operator of E, denoted 
by £\ is that linear partial differential operator mapping rt-index tensor fields into m-index tensor 
fields such that 

g'^-m.^.^u... - [ft(gp--)]^--/M.... = v^J^ (1) 

where is some vector field depending on the fields / and g. From Eq. (1) we can see that this 
definition automatically guarantees that, if the field / is a solution of the linear system £{J) — and 
g a solution of the adjoint system ft (g) ~ 0, then we obtain the local continuity law Vf^J'^ = 0, which 
establishes that J'^ is a covariantly conserved current. This fact means that for any homogeneous 
equation system, one can always construct a conserved current taking into account the adjoint 
system. This general result contains the self-adjoint case {£ = £'') as a particular one, for which / 
and g correspond to two independent solutions. 

In the present work, / and g will be associated with the first-order variations of the backgrounds 
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fields. Such field variations will correspond, on the phase space, to one-forms [1]. In this manner, 
the left-hand side of Eq. (1) can be understood as a wedge product on such phase space: g A £(/) — 
^'^{9) ^ / = ^tJ.J^, and something similar for the bilinear form in its dependence on the fields / 
and g (the operators £, £\ and are depending only on the background fields). 

Traditional branch: It is important, for completeness, that we outline the original idea for 
introducing the definition (1) in Ref. |Q: reduction of systems of linear partial differential equations 
to equations for scalar potentials (called Debye potentials), which determine a complete solution of 
the original system. 

If we have the linear system £{f) = 0, and there exist linear operators such that 

S£ = OT, 

identically, then the field S''{ip) satisfies the equation 

£t(5t(V))=0, 

provided that the scalar field satisfies 

In particular, if £ is self-adjoint {£^ ~ £), then / = S^{^/j) is a solution of £{f) = 0. For example, 
in the case considered in the present work, the (matrix) operator governing the field perturbations 
in the Einstein-Maxwell and Einstein- Weyl theories are, in fact, self-adjoint ^. 

Moreover, the existence of operators S, O, and T satisfying the above identity, is equivalent to 
the existence of a decoupled system 

0(*) = 0, 

obtained from the original system £{f) — 0, such that the scalar field ^' = T{f). 

Now, we can mix both branches of the adjoint operators scheme: since the fields V' and 5* satisfy 
equation which are adjoints to each other, we can establish, in according to the first branch, that 

^-(0*) - (oV)* - v^j^(Va*), 

which means that J'^('0, ^') = 0. Furthermore, since *is finally depending on (^f = r(/) = r(5t(-0))), 
is dependent only on -0 (however, see the paragraph before Eq. (8)). 
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On the other hand, ahhough this result has been estabhshed assuming only tensor fields and the 
presence of a single equation, such a result can be extended in a direct way to equations involving 
spinor fields, matrix fields, and the presence of more than one field 

It is worth pointing out some issues on the notation. The first-order field variations appearing 
in Ref. [5,6] are denoted by a superscript B. On the other hand, the field variations coincide, in 
according to Witten's interpretation [0, with an infinite-dimensional generalization of the usual 
exterior derivative, which is traditionally represented by the symbol 6. However, in Ref. [5,6], the 
Newman-Penrose formalism is used, in which the symbol S is employed for denoting one of the direc- 
tional derivatives defined by the null tetrad. In this manner, for avoiding confusion, we will maintain 
the symbol S as usual in the Newman-Penrose notation, and the superscript B for the first-order 
field variations (the exterior derivative of background fields). In the present article, the exterior 
derivative will not be performed explicity, and it will be sufficient for our purposes to understand 
any quantity with the superscript B as a one-form on the phase space. Quantities without such a su- 
perscript will correspond to background fields, which mean zero- forms on the phase space Q]. With 
these previous considerations, formulae and notation of Ref. [5,6] will be used throughout this paper. 



In Ref. [g, the metric and the vector potential variations of the Einstein-Maxwell equations 
have been obtained starting from a decoupled set of equations governing certain combinations of 
field variations, and have the form: 



where the 4x4 matrix operator O can be read from Eq. (24) of that reference, and (^'^) is the 



III. COVARIANT CONTINUITY LAWS 



Oi^'^) = 0, 



(2) 



matrix: 



(3) 
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where 

= (2(/.i^'f -3^'2</)^)/2,/)i, 

= n^ + {D-2e-p){^f/2<l>,), (4) 

= a'' + iS-2p-T){cj)^/2(f>,). 

The decoupled system (2) was obtained provided that the background electromagnetic field is non- 
null with one of its principal null directions being geodetic and shear free, which means that ^'o = 
= ^*1, K ^ — a, and 00 = in the background Note that the variations of these vanishing 
background quantities are those involved in ('5^). As described in previous Section, the variations of 
background fields correspond to one-forms, in particular, those ones in Eq. (4). However, as shown 
in Appendix B, they correspond to exact one-forms. This property will be useful below. 

On the other hand, according to Wald's method used in Ref. the matrix potential (V'), which 
determines the complete variations, satisfies the equations (see Eqs. (25) and (27) of Ref. [^) 

O^(^)=0, (5) 

where 



(6) 



/ Ml, ^^ 

-Mo 

V -V'E / 

Thus, from the definition (1) we have that (i/') A 0{^^) ~ 0^'^P) A {^^) = J^(*^, i/-). The left- 
hand side contains terms of the form My A Oii(*o ) - C'|i(Afi/) A , where On = S-4:a + Tr and 
oil = ((5 — 4q; + 7r)t = S + 3a + f3 — T, are the entries 11 of the 4x4 matrices O and O'' respectively 
(see Eqs. (24) and (25) in Ref. Q). Taking into account that 6 = m^9^, and that is acting on scalar 
fields, such a term can be expressed in the form My AOn{'i'o)~0\i{Mv)A'i'^ = V^(m^Mi' A^'Jf), 
and similarly for the remaining terms considering that the tetrad vectors D, A, and S are defined 
by I'^dp,, n^df^, and m^df^ respectively. In this manner, from Eqs. (2) and (5) we obtain cxplicity 
the following covariant conservation law: 

= 0, Jf = -li'iMy A + 201 a a^) - n^'iMo' A ^Jf + 20iVe A k^) 

+m^'{Mo' A *f + 20ii/'g A k^) + to^(Mi, A + ^iV-e A ct^). (7) 
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It is very easy to determine that, like (^'^) in Eq. (3), in Eq. (6) is made out of one-forms. 
From Eqs. (28) and (29) of Ref. j^, the metric and vector potential variations are given by 

= 5t(V'), 

where 5^ is an operator depending only on background fields (see Eqs. (26), (28), and (29) in Ref. 
1^), which correspond to zero-forms. Since /i^^ and correspond to one-forms, from the preceding 
equation, (ip) correspond also to one-forms. This implies automatically that in Eq. (7) is a (non- 
degenerate) two-form on the corresponding phase space. However, we have to demonstrate that 
corresponds to a closed two-form for obtaining a symplectic structure on the phase space. For this 
purpose, we use the results in Appendix B. For example, the first term — Z^A/i' A '^f for J'' in Eq. 
(7) can be rewritten as: 

where Eq. (Bl) has been used; the above expression means that such a term is an exact two- form 
(and automatically a closed two-form). Similarly, using Eqs. (Bl) and (B2), the remaining terms 
in Eq. (7) correspond to closed two-forms. Therefore, itself is a closed two-form, as required for 
constructing a symplectic structure. Such geometrical structure is defined hy u = J-^ J^dE^, where 
S is an initial value hypersurface. As J'^ is conserved, lu is independent of the choice of T, and, in 
particular, is Poincare invariant 0. In this manner, cj is a covariant symplectic structure on the 
phase space. 

On the other hand, as shown in the Appendix A, the one- forms 'i>^, k^, and a^, are 
invariant under gauge transformations of the vector potential variations, h^^ ^ + V^Xj ^.nd thus, 
J'' and Lo have the same invariance property (in contrast with those of Refs. [1-4]). In field space, 
V^X corresponds to the gauge directions, therefore, this invariance property defines automatically 
a gauge- invariant two- form uj on the subtler space Z = Z/G, being Z the space of solutions and 
G the group of gauge transformations [0. Similarly, since such one- forms are also invariant under 
infinitesimal rotations of the tetrad (see Appendix A), we have also a two- form lu independent on 
the perturbed tetrad gauge freedom, defined on the subtler space T = Z/H, being H the group of 
infinitesimal rotations of the perturbed tetrad. However, J'^ and u are not gauge invariant with 
respect to infinitesimal diffeomorphisms (in similarity to those of Refs. [1-4]). 
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From Eq. (7) we can see that the conserved current depends on the background fields and the 
solutions admitted by the decoupled system and its adjoint. However, the entries of the matrix (^'^) 
can be defined in terms of the metric and vector potential variations, which in turn, are defined in 
terms of the potential matrix {tp) (see Eqs. (28) and (29) of Ref. |^), and therefore, and lu can be 
expressed in terms of a single complex solution of the equations for the potential (ip). However, in 
the more general case, if (ip)! is a solution admitted by the equations for the potentials, the matrix 
(^f^) can be expressed in terms of a second solution {tp)2, in general different of (?/')i- Therefore, 

and Lo are expressed in terms of a pair of solutions for those equations. Finally, J'^ and lu will be 
given completely in terms of the potentials (one- forms), which become the fundamental geometrical 
structures. In this manner, the analysis of the structure of the phase space has been reduced to the 
study of equations for the potentials, which is a relatively simple issue. 

Similarly, when the background electromagnetic field is null (and there is a possibly nonzero 
cosmological constant) the continuity equation for the field variations in the Einstein-Maxwell theory 
is given by (see Eqs. (21) and (24) of Ref. §) 

V^,r^O, .r = l'' [V'G A ( A - 47 + Ai)*J^ + A ( A - 27 + fi)^^ 
+2^2 (V^G A{D- 2e)(^^ -{D + ie + 3p)yjG A ^^)] 
-n^[(D + 4e + 3p)Vg A + (Z) + 2e + p)V'E A iff] 
-m^[7/;G A (^ - 4a + Tr)^'^ + tpE A (S ~ 2a + n)(p^] 
-t-m^ [(5 + 4/3 + 3r)V'G A + ((5 + 2/3 + r)V'E A (p^] , (8) 

where similarly is a closed two-form on the phase space, and the corresponding symplectic 
structure uj = J^dS^ is Poncaire invariant because is conserved. Furthermore, and the 
corresponding symplectic form uj have essentially the same invariance properties that the previous 
case (see Appendix A), and therefore are defined in a gauge-invariant form on the corresponding 
subtler spaces Z and T. Analogously, they are given, in general, in terms of a pair of solutions 
admitted by the corresponding equations for the potentials (see Eqs. (24) in Ref. [^). 

In the case of solutions of the Einstein- Weyl equations such that the flux vector of the neutrino 
field is tangent to a shear-free congruence of null geodesies, the field variations satisfy the continuity 
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equation (see Eqs. (38)-(47) of Rcf. [g): 

Vf.ll'^i-Mi- A ^-f - rjitpG A (T^ - ik-qiriyMi> A + 2ikr]ir]yMo' A 5-^] 

-n^[Mo' A^-^ + r/i^N Ak^] 
+m^[r;iV'G A + Mo' A - ifcr^i^^Afo' A k^] + m^[Mi. A + mV'N A ct^]} = 0. (9) 

where Mq', My, ip]^, and ipc are the potentials, and similarly \I>q, Xf^^, ct^, and are defined in 
terms of such potentials (see Eqs. (47) of Ref. |^). 

The continuity equations (7)-(9) deserve some additional comments. The existence of such 
continuity laws follows directly from the definition of the adjoint of a differential operator, and we 
have only to identify conveniently the system of differential equations for applying such definition. 
In this sense, the adjoint operators scheme allows us to reveal the underlying geometrical structures 
for the corresponding theory, in a elementary and straightforward way. The argument involved is, 
as seen above, very simple. 

On the other hand, in the Lagrangian approach known in the literature ||^ , the symplectic cur- 
rents arise from the exterior derivative of the so called symplectic potentials, which correspond to 
one-forms on the phase space. This guarantees that such bilinear forms correspond automatically to 
closed two-forms (symplectic structures) on the phase space. However, in the present approach, the 
symplectic currents in Eqs. (7)"(9) are constructed in a direct way using only the concept of adjoint 
operators, regardless of the existence of symplectic potentials. Furthermore, our symplectic currents 
are generated from the usual derivatives (essentially directional derivatives) of the potentials, which 
correspond also to one-forms on the phase space. This suggests a possible relationship between the 
fundamental structures of the different approaches: the symplectic potentials of the Lagrangian ap- 
proach and the potentials considered in the present scheme. This may be the subject of forthcoming 
works. 



V. CONCLUDING REMARKS 



As we have seen, the present approach based on the concept of the adjoint of a differential operator. 
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represents a more direct and elementary procedure for constructing geometrical structures manifestly 
covariant on the phase space. Specifically, it has allowed us to demonstrate that, for the restricted 
class of spacetimes considered here, there exist two covariantly conserved currents for the field 
variations of each theory, the symplectic current (which coincides fully with that obtained from the 
Lagrangian approach §), and that found in the present approach (which is no obtained directly 
from a Lagrangian approach). Since both currents can be expressed, finally, in terms of solutions 
of the equations for the corresponding scalar potentials, our fundamental geometrical structures, a 
direct comparison is possible and convenient in order to establish the differences and similarities. 
On the other hand, although some restrictions on the background fields have been required, a lot 
of solutions of wide interest is contained in the cases considered here: rotating (charged) black 
holes, colliding plane waves, etc. The detailed study of the corresponding currents and geometrical 
structures associated with those particular solutions will be the subject of subsequent works. In 
fact, as mentioned in the Introduction, works along these lines have permitted to demonstrate 
that, the continuity equation (8) directly yields, in the particular case of the gravitational and 
electromagnetic perturbations of the Schwarzschild black hole, the conservation of the energy for 
those perturbations without invoking, as usually, the Schrodinger-type equation for establishing 
such a conservation relation. This result is connected with the existence of conserved quantities for 
non-Hcrmitian system of differential equations . 

Finally, the general result established in Sec. II can be understood as an important extension of 
the original Wald's method wherever there exists an appropriate decoupled equation, it is not 
only possible to express the complete solution in terms of scalar potentials, but also to find auto- 
matically a corresponding (covariantly) conserved current, which can be expressed, in turn, in terms 
of the same potentials. The possible applications of such a result in the modern theories involv- 
ing gravity and matter fields such as string theory (and other areas of physics) are open questions . 
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APPENDIX A: INVARIANCE PROPERTIES 

Although it is not mentioned explicitly in Refs. ^, the perturbed quantities appearing in the 
decoupled equations are invariant under infinitesimal rotations of the tetrad. For example, in the 
case of the perturbed quantity (see Eqs. (4)) made out of the perturbations of the spin coefficient 
K and of the component of electromagnetic field 0Oj we have that: 

= m'^m^{Tl,r -]^{^-T)m-h^,-{D-2e-pWml), (Al) 

where (F^^)^ — \g'^''{V ^h^p + V ^h^p — S/ ph^^) and /i^^ is the metric perturbation. In this manner, 
the first and second terms in the above equation are defined completely in terms of hpy. On the 
other hand, l^rn^ is dependent on the perturbed tetrad gauge freedom. Furthermore, from the 
definition (f)Q = l^m'^Fp„, one finds that 

^ _ +^)/^^m''V + 20i(?^m;^), (A2) 

where 0i is the only nonvanishing component of the background electromagnetic field, and F^^^ — 
dpbi, — d^bp, being bp the vector potential perturbations. Therefore, from Eqs. (Al) and (A2) we can 
see easily that the quantity — + {D ~ 2e — p)((/)q /20i) is independent of the perturbed tetrad 
gauge freedom. Similarly for the remainding quantities ^'q,^'^, and in Eq. (4). Thus, (^'^) is 
invariant under the transformation considered. Since the field perturbation Fj^^ is invariant under 
the ordinary gauge transformation bp bp+VpX, where x is an arbitrary scalar field, 0q in Eq. (A2) 
is also invariant under this transformation and, in this manner, (^'^) is (see Eqs. (4)). However, 
it is not difficult to demonstrate that (^'^) is not invariant under infinitesimal diffeomorphisms 
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bfj. ^ bfj_ + C^A^, and /i^i^ ft,^^ + C^g^v, where is an arbitrary vector field, and and gf^v 
the background fields. 

Similarly, in the case of solutions of the Einstein-Maxwell theory with a null background electro- 
magnetic field {4)2 is the only nonvanishing component) Q, we find that 

00 = I'^TT^'F]^. + \¥2l^m''h^,, (A3) 

which is independent of the perturbed tetrad gauge freedom and, of course invariant under — > 
6^ + VpX; however it is not under infinitesimal diffeomorphisms. On the other hand, it is easy to 
verify that 'I'g is invariant under ordinary gauge transformations and infinitesimal diffeomorphisms 
(and independent of the perturbed tetrad gauge freedom). In this manner, our quantities are, in 
general, only invariant under both ordinary gauge transformations of the vector potential pertur- 
bartions and rotations of the tetrad. 



APPENDIX B: EXACT ONE-FORMS ON THE PHASE SPACE 



The field variations appearing in the decoupled equation (3) correspond to exact one-forms on 
the phase space. 

The case of trivial: \l/g is directly the exterior derivative of the zero-form \l/o (see last 

paragraph in Sec. II). 

The other cases are more elaborated. For example, VE"^ can be rewritten as: 



since (f)Q = Q in the background. Similarly, k , and a can be rewritten as: 



a = a 



201 _ 



(B2) 
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Therefore, from Eqs. (Bl) and (B2), "iif, k^, and correspond to the exterior derivative of 
background zero-forms, and then correspond to exact one-forms. A consequence of this property is 
that (^'^) is made out of closed one- forms: 

= 0. (B3) 

Similarly, the field variations ^'q , ^t^, ii^, and appearing in Eq. (9) (and defined in Eqs. (39) 
of Ref. dl), correspond to exact one- forms on the phase space. 

The cases of and 0q in Eq. (8) are also trivials: they are directly the variations of the corre- 
sponding background quantity. 
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